arXiv: 1505.07812v2 [nucl-ex] 7 Nov 2015 


Azimuthal anisotropy in U+U and An+An collisions at RHIC 


L. Adamczyk/ J. K. Adkins,G. Agakishiev,^® M. M. Aggarwal,®^ Z. Ahammed,^^ 1. Alekseev,^® J. Alford,^® 

A. Aparin,^® D. Arkhipkin,® E. C. Aschenauer,® G. S. Averichev,®® A. Banerjee,^®’ R. Bellwied,^® A. Bhasin,®'’ 

A. K. Bhati,®® P. Bhattarai,^^ J. Bielcik,®® J. Bielcikova,®® L. G. Bland,® I. G. Bordyuzhin,®® J. Bouchet,®® 

A. V. Brandin,^® 1. Bunzarov,®® T. P. Burton,® J. Butterworth,®® H. Caines,®® M. Calderon de la Barca Sanchez,® 
J. M. Campbell,^® D. Cebra,® M. C. Cervantes,^® I. Chakaberia,® P. Chaloupka,®® Z. Chang,^® S. Chattopadhyay,"®®’ 

J. H. Chen,®® X. Chen,®® J. Cheng,'®^ M. Cherney,® W. Christie,® G. Contin,®® H. J. Crawford,^ S. Das,®® 

L. C. De Silva,® R. R. Debbe,® T. G. Dedovich,®® J. Deng,®® A. A. Derevschikov,®® B. di Ruzza,® L. Didenko,® 

C. Dilks,®® X. Dong,®® J. L. Drachenberg,®® J. E. Draper,® C. M. Du,®® L. E. Dunkelberger,® J. C. Dunlop,® 

L. G. Efimov,®® J. Engelage,^ G. Eppley,®® R. Esha,® O. Evdokimov,® O. Eyser,® R. Fateini,®® S. Fazio,® 

P. Federic,®® J. Fedorisin,®® Z. Feng,'’ P. Filip,®® Y. Fisyak,® C. E. Flores,® L. Fulek,® C. A. Gagliardi,^® 

D. Garand,®® F. Geurts,®® A. Gibson,'®® M. Girard,^® L. Greiner,®® D. Grosnick,'®® D. S. Gunarathne,'®® 

Y. Guo,®®" S. Gupta,®®' A. Gupta,®® W. Guryn,® A. Hamad,®®* A. Hamed,'®® R. Haque,®®' J. W. Harris,®® L. He,®® 

S. Heppelmann,® S. Heppelmann,®® A. Hirsch,®® G. W. Hoffmann,^® D. J. Hofman,® S. Horvat,®® H. Z. Huang,® 

B. Huang,® X. Huang,^4 P. Huck,®" T. J. Humanic,®® G. Igo,® W. W. Jacobs,®® H. Jang,®® K. Jiang, ®®' E. G. Judd,4 
S. Kabana,®® D. Kalinkin,®® K. Kang,44 K. Kauder,4® H. W. Ke,® D. Keane,®® A. Kechechyan,®® Z. H. Khan,® 

D. P. Kikola,4® L Kisel,®® A. Kisiel,^® D. D. Koetke,4® T. Kollegger,®® L. K. Kosarzewski,4® L. Kotchenda,®® 

A. F. Kraishan,4® P. Kravtsov,®® K. Krueger,® L Kulakov,®® L. Kumar,®® R. A. Kycia,®® M. A. C. Lamont,® 

J. M. Landgraf,® K. D. Landry,® J. Lauret,® A. Lebedev,® R. Lednicky,®® J. H. Lee,® W. Li,®® Y. Li,44 C. Li,®®” 

Z. M. LiJ X. Li,4® X. Li,® M. A. Lisa,®® F. Liu,®^ T. Ljubicic,® W. J. Llope,4® M. Lomnitz,®® R. S. Longacre,® 

X. Luo,®' L. Ma,®® R. Ma,® Y. G. Ma,®® G. L. Ma,®® N. Magdy,®° R. Majka,®® A. Manion,®® S. Margetis,®® 

C. Markert,4® H. Masui,®® H. S. Matis,®® D. McDonald,4® K. Meehan,® N. G. Minaev,®® S. Mioduszewski,4® 

B. Mohanty,®®' M. M. Mondal,4® D. A. Morozov,®® M. K. Mustafa,®® B. K. Nandi, ®4 Md. Nasim,® T. K. Nayak,4®' 

G. Nigmatkulov,®® L. V. Nogach,®® S. Y. Noh,®® J. Novak,®® S. B. Nurushev,®® G. Odyniec,®® A. Ogawa,® K. Oh,®4 
V. Okorokov,®® D. L. Olvitt Jr.,40 B. S. Page,® R. Pak,® Y. X. Pan,® Y. Pandit,® Y. Panebratsev,®® B. Pawlik,®® 

H. Pei,®" C. Perkins,4 A. Peterson,®® P. Pile,® M. Planinic,®® J. Pluta,4® N. Poljak,®® K. Poniatowska,4® J. Porter,®® 

M. Posik,4® A. M. Poskanzer,®® N. K. Pruthi,®® J. Putschke,4® H. Qiu,®® A. Quintero,®® S. Ramachandran,®® 

S. Raniwala,®® R. Raniwala,®® R. L. Ray,®® H. G. Ritter,®® J. B. Roberts,®® O. V. Rogachevskiy,®® J. L. Romero,® 
A. Roy,4®' L. Ruan,® J. Rusnak,®® O. Rusnakova,®® N. R. Sahoo,®® P. K. Sahu,®® 1. Sakrejda,®® S. Salur,®® 

J. Sandweiss,®® A. Sarkar,®® J. Schambach,®® R. P. Scharenberg,®® A. M. Schmah,®® W. B. Schmidke,® 

N. Schmitz,®® J. Seger,® P. Seyboth,®® N. Shah,® E. Shahaliev,®® P. V. Shanmuganathan,®® M. Shao,®®' B. Sharma,®® 

M. K. Sharma,®® W. Q. Shen,®® S. S. Shi,® Q. Y. Shou,®® E. P. Sichtermann,®® R. Sikora,® M. Simko,®® 

M. J. Skoby,®® D. Smirnov,® N. Smirnov,®® L. Song,®® P. Sorensen,® H. M. Spinka,® B. Srivastava,®® 

T. D. S. Stanislaus,®® M. Stepanov,®® R. Stock,®® M. Strikhanov,®® B. Stringfellow,®® M. Sumbera,®® B. J. Summa,®® 
X. Sun,®® X. M. Sun,7 Z. Sun,®® Y. Sun,®^ B. Surrow,®® D. N. Svirida,®® M. A. Szelezniak,®® Z. Tang,®^ 

A. H. Tang,® T. Tarnowsky,®® A. N. Tawfik,®® J. H. Thomas,®® A. R. Timmins,®® D. Tlusty,®® M. Tokarev,®® 

S. Trentalange,® R. E. Tribble,®® P. Tribedy,®®' S. K. Tripathy,®® B. A. Trzeciak,®® O. D. Tsai,® T. Ullrich,® 

D. G. Underwood,® 1. Upsal,®® G. Van Buren,® G. van Nieuwenhuizen,® M. Vandenbroucke,®® R. Varma,®® 

A. N. Vasiliev,®® R. Vertesi,®® F. Videbaek,® Y. P. Viyogi,®®" S. Vokal,®® S. A. Voloshin,®® A. Vossen,®® F. Wang,®® 
Y. Wang,®® H. Wang,® J. S. Wang,®® Y. Wang,®" G. Wang,® G. Webb,® J. C. Webb,® L. Wen,® G. D. Westfall,®® 

H. Wieman,®® S. W. Wissink,®® R. Witt,®® Y. F. Wu,®’ Z. Xiao,®® W. Xie,®® K. Xin,®® Y. F. Xu,®® N. Xu,®® Z. Xu,® 
Q. H. Xu,®® H. Xu,®® Y. Yang,®’ Y. Yang,®® C. Yang,®®’ S. Yang,®®' Q. Yang,®®’ Z. Ye,® P. Yepes,®® L. Yi,®® K. Yip,® 

I. -K. Yoo,®® N. Yu,®’ H. Zbroszczyk,®® W. Zha,®®’ X. P. Zhang,®® J. B. Zhang,®’ J. Zhang,®® Z. Zhang,®^ S. Zhang,®^ 
Y. Zhang,®®’ J. L. Zhang,®® F. Zhao,® J. Zhao,®^ C. Zhong,®® L. Zhou,®®^ X. Zhu,®® Y. Zoulkarneeva,®® and M. Zyzak®® 

(STAR Collaboration) 

^AGH University of Science and Technology, Cracow 30-059, Poland 
^Argonne National Laboratory, Argonne, Illinois 60439, USA 
^ Brookhaven National Laboratory, Upton, New York 11973, USA 
‘^University of California, Berkeley, California 94720, USA 
^University of California, Davis, California 95616, USA 
® University of California, Los Angeles, California 90095, USA 
^Central China Normal University (HZNU), Wuhan 430079, China 


2 


® University of Illinois at Chicago, Chicago, Illinois 60607, USA 
^Creighton University, Omaha, Nebraska 68178, USA 
Czech Technical University in Prague, FNSPE, Prague, 115 19, Czech Republic 
Nuclear Physics Institute AS CR, 250 68 Rez/Prague, Czech Republic 
Frankfurt Institute for Advanced Studies FI AS, Frankfurt 60438, Germany 
Institute of Physics, Bhubaneswar 751005, India 
^'^Indian Institute of Technology, Mumbai 400076, India 
Indiana University, Bloomington, Indiana 4'^408, USA 
Alikhanov Institute for Theoretical and Experimental Physics, Moscow 117218, Russia 
University of Jammu, Jammu 180001, India 
Joint Institute for Nuclear Research, Dubna, Ifl 980, Russia 
Kent State University, Kent, Ohio 44^4^! USA 
University of Kentucky, Lexington, Kentucky, 40506-0055, USA 
Korea Institute of Science and Technology Information, Daejeon 305-701, Korea 
Institute of Modem Physics, Lanzhou 730000, China 
Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
Max-Planck-Institut fur Physik, Munich 80805, Germany 
Michigan State University, East Lansing, Michigan 48824, USA 
Moscow Engineering Physics Institute, Moscow 115409, Russia 
National Institute of Science Education and Research, Bhubaneswar 751005, India 
Ohio State University, Columbus, Ohio 43210, USA 
Institute of Nuclear Physics PAN, Cracow 31-342, Poland 
Pan jab University, Chandigarh 160014, India 
^^Pennsylvania State University, University Park, Pennsylvania 16802, USA 
Institute of High Energy Physics, Protvino 142281, Russia 
^^Purdue University, West Lafayette, Indiana 4^907, USA 
Pusan National University, Pusan 609735, Republic of Korea 
University of Rajasthan, Jaipur 302004, India 
Rice University, Houston, Texas 77251, USA 
University of Science and Technology of China, Hefei 230026, China 
Shandong University, Jinan, Shandong 250100, China 
Shanghai Institute of Applied Physics, Shanghai 201800, China 
Temple University, Philadelphia, Pennsylvania 19122, USA 
Texas A&M University, College Station, Texas 77843, USA 
University of Texas, Austin, Texas 78712, USA 
University of Houston, Houston, Texas 77204, USA 
‘^'^Tsinghua University, Beijing 100084, China 
United States Naval Academy, Annapolis, Maryland, 21402, USA 
Valparaiso University, Valparaiso, Indiana 46383, USA 
Variable Energy Cyclotron Centre, Kolkata 700064, India 
Warsaw University of Technology, Warsaw 00-661, Poland 
Wayne State University, Detroit, Michigan 48201, USA 
World Laboratory for Cosmology and Particle Physics (WLCAPP), Cairo 11571, Egypt 
Yale University, New Haven, Connecticut 06520, USA 
University of Zagreb, Zagreb, HR-10002, Croatia 

Collisions between prolate uranium nuclei are used to study how particle production and azimuthal 
anisotropies depend on initial geometry in heavy-ion collisions. We report the two- and four- particle 
cumulants, U2{2} and U2{4}, for charged hadrons from U-l-U collisions at ^snn = 193 GeV and 
Au-|-Au collisions at ^snn = 200 GeV. Nearly fully overlapping collisions are selected based on 
the energy deposited by spectators in Zero Degree Calorimeters (ZDCs). Within this sample, the 
observed dependence of V2{2} on multiplicity demonstrates that ZDC information combined with 
multiplicity can preferentially select different overlap configurations in U-l-U collisions. We also 
show that V 2 vs multiplicity can be better described by models, such as gluon saturation or quark 
participant models, that eliminate the dependence of the multiplicity on the number of binary 
nucleon-nucleon collisions. 

PACS numbers: 25.75.-q, 25.75.Ag, 25.75.Ld 

Collisions of nuclei at the Relativistic Heavy-Ion Col¬ 
lider (RHIC) and the Large Hadron Collider (LHC) cre¬ 
ate a fireball hot and dense enough to form a Quark 
Gluon Plasma (QGP) [H- Anisotropies in the final mo¬ 


mentum space distributions can be traced back to spatial 
anisotropies in the initial state and are used to under¬ 
stand the nature of the fireball These anisotropies 

are studied using harmonics of the distribution of the az- 
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imuthal angle (f) separation between pairs of particles [3“ 
1^. The inference of the properties of the fireball from 
these measurements is limited however by uncertainties 
in the description of the initial state Q- Collisions be¬ 
tween uranium nuclei, which have an intrinsic prolate 
shape Q, provide a way to manipulate this initial ge¬ 
ometry to test our understanding of the initial state of 
heavy-ion collisions and the subsequent fireball Q. 

Even in nearly fully overlapping collisions of U nuclei 
(impact parameter 6 Ri 0 fm), the initial matter distribu¬ 
tion can exhibit very different shapes. In one extreme, 
the major axes of both colliding nuclei could lie parallel 
to the beam so that the tip of one nucleus impinges on 
the tip of the other (tip-tip). Another extreme occurs if 
the major axes of the nuclei are parallel to each other but 
perpendicular to the beam so that they collide side-on- 
side or body-body. There are two principal differences in 
these two configurations — tip-tip collisions have a larger 
number of binary nucleon-nucleon collisions A^bin while 
body-body collisions have a smaller Abin but a more el¬ 
liptic overlap region (larger eccentricity £2). The larger 
Abin in the tip-tip configuration is expected to lead to 
a larger multiplicity of produced particles [l^ . 11 | while 
the more elliptic shape of the body-body collisions is ex¬ 
pected to lead to a larger second harmonic anisotropy V 2 ■ 
The dependence of V 2 on multiplicity in nearly fully over¬ 
lapping U-l-U collisions therefore tests our understand¬ 
ing of particle production and the development of V 2 - 
An anti-correlation between V 2 and multiplicity in these 
collisions will also demonstrate that multiplicity can be 
used to select enhanced samples of body-body or tip- 
tip configurations. Those samples can then be used to 
study other topics like the path-length dependence of 
jet-quenching [9|, or the extent to which three-particle 
charge-dependent correlations BP can be attributed 
to local parity violation [l^ or background effects 0. 
We also investigate two models that do not include any 
explicit dependence on A^bin; one based on gluon satu¬ 
ration B, Bl and the other based on the number of 
participating constituent quarks B0- 

In this Letter we report measurements of the two- 
and four-particle cumulant of V 2 {v 2 { 2 } and f2{4}) in 
197 Au-h197Au and -b collisions at = 200 

and 193 GeV respectively. Both minimum bias and 
nearly fully overlapping events where most of the nu¬ 
cleons participate in the collision are studied. The data 
sets were collected by STAR Bl ia 2011 and 2012. The 
U-bU data consists of approximately 307 million events 
including 7 million specially triggered central events. 
Charged particles within pseudo-rapidity window |?7| < 1 
were detected using the STAR Time Projection Cham¬ 
ber (TPC) Bl- We select tracks within the transverse 
momentum range 0.2 < pr < 2.0 GeV/c. The STAR 
Zero-degree Calorimeters (ZDCs) B were used to select 
the sample of nearly fully overlapping events; those hav¬ 
ing large multiplicity but little activity in the ZDCs. The 


ZDC resolution was determined to be 23 ± 2% from the 
observation of the single neutron peak in the ADC signal. 
The ZDC selection requires ZDCs on both sides of the 
detector to have a signal smaller than the specified cut. 
The tracking efficiency is corrected via embedding and 
weights in 77 and (p derived from the inverse of the distri¬ 
bution of tracks observed over many events. This method 
allows us to correct our 772{2,4} measurements for imper¬ 
fections in the tracking efficiency 712{4} was calculated 
using the Q-Cumulant method [23| while 7;2{2} was cal¬ 
culated directly from particle pairs (cos2((;iii — 4'2))- To 
reduce the contribution from HBT, Coulomb and track¬ 
merging effects, a minimum rj separation of jAryj > 0.1 
is required for 7i2{2}. Measurement uncertainties were 
estimated by varying event and track selection criteria, 
varying efficiency estimates, and by comparing data from 
different run periods. These uncertainties are quite small; 
less than 0.1% absolute variation on 7;2{2,4}. 

Figure [T] shows the two- and four- particle cumulant 
7;2{2} and 7'2{4} from minimum bias 200 GeV Au-bAu 
and 193 GeV U-bU collisions as a function of efficiency 
corrected charged particle multiplicity dNch/di]. We 
find that the relationship between dNch/dr] and central¬ 
ity fraction can be parameterized as {dNch/drjy^'^ = 
Cl — C 2 X + C3 exp(— 040;“'^) with Ci=5.3473, C2=4.298, 
C3=0.2959, C4=18.21, and cg = 0.4541 for U-bU and 
ci=5.0670, C2=3.923, C3=0.2310, C4 = 18.37, and 
C5=0.4842 for Au-bAu. Multiplicity trends for 7’2{2} and 
772 {4} in U-bU collisions are mostly similar to those ob¬ 
served in Au-bAu collisions. A notable difference however 
is seen in the 7>2{4} measurements in central U-bU col¬ 
lisions. Whereas 7;2{4} (shown in the inset) is negative 
for central Au-bAu collisions, it is positive for U-bU col¬ 
lisions. Previous studies showed that fluctuations in the 
number of participating nucleons cause 74 {4} in central 
Au-bAu collisions to become negative [2^. The obser¬ 
vation of 7;|{4} > 0 in the most central U-bU collisions 
indicates that the prolate shape of uranium increases the 
anisotropy in the final momentum space distributions of 
the observed particles. 

Glauber-based models have typically used a two- 
component model ((l-a;hard)A^part/2-bXhardA^bin) for the 
multiplicity, where Apart is the number of struck nucle¬ 
ons, Abin is the number of binary nucleon-nucleon colli¬ 
sions, and cchard is a fractional contribution of Abin to the 
multiplicity B; Bl • The multiplicity is then assumed to 
fluctuate according to a convolution of negative binomial 
distributions (NBD) with parameters n and k related to 
the mean and width measured from p+p collisions at the 
same energy and in the same \ri\ window (B - We will re¬ 
fer to this model as “Glauber-Xhard-” Since the number of 
hard scatterings is known to scale with Abin, a^hard is of¬ 
ten thought of as reflecting the contribution of hard pro¬ 
cesses to the multiplicity. It can also be thought of as a 
coherence parameter with Xhard = 1 giving the maximum 
incoherence as multiplicity entirely arises from indepen- 
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FIG. 1. (Color online) The two- and four- particle cumulant 
ii 2 { 2 } and U 2 { 4 } within I 77 I < 1 versus dNch/dri from 200 GeV 
Au-|-Au and 193 GeV U-l-U collisions. Dashed lines show 
U-fU centralities based on dNch/drj measured in I 77 I < 0.5. 
u|{4} (the experimentally observed quantity) is shown in the 
inset without taking the fourth root in the range where it is 
near zero or negative. 


dent binary nucleon-nucleon collisions. The Glauber- 
a^hard model indicates that V 2 in U-l-U collisions should 
begin to decrease markedly for events with multiplicities 
in the top 1% [13| forming a knee structure where tip- 
tip collisions with larger iVbin and smaller eccentricity 
begin to dominate. Vertical dashed lines in the figure in¬ 
dicate the 1%, 0.1% and 0.01% highest multiplicity U-l-U 
collisions. No knee structure is observed suggesting the 
Glauber-cchard model may not be the correct description. 
Adding more multiplicity fluctuations causes the knee 
structure to disappear 26| but this will also significantly 
increase the average £2 in central collisions. 

To explore the dependence of V 2 on the initial eccen¬ 
tricity £ 2 , we plot V 2 /e 2 versus dNch/drj. It was found 
previously that V 2 IS 2 monotonically increases with in¬ 
creasing dNchjdr] and depending on the model for the 
initial eccentricity may, or may not saturate in the most 
central collisions [ 2 ^. Figure [2] shows U2{2}/£2{2} and 
'y2{4}/£2{4} from Au-bAu and U-l-U collisions. £2{2} 
and £2{4} are the second and fourth cumulants of the 
participant eccentricity distributions calculated from the 
Glauber-cchard model [27l - E9j| . Both U-l-U and Au-|-Au 
follow a similar trend for t' 2 /£ 2 - However, a turn-over is 
observed in central collisions {dNch/drj > 500). This has 
not been observed previously since measurements have 
typically been integrated over 5% most central [^. The 
turn-over is consistent with the model overestimating £2 
in central collisions. Increasing the multiplicity fluctu¬ 
ations as in Ref. [ 2 ^ will only increase the eccentricity 
in central collisions suggesting that a different explana¬ 
tion may be required to explain both the turn-over of 
V 2 IS 2 and the lack of a knee structure in V 2 vs dNch/drj. 
Using a new set of Woods-Saxon parameters derived 


FIG. 2. (Color online) V 2 scaled by participant eccentricity 
from 200 GeV Au-|-Au and 193 GeV U-fU collisions. The 
eccentricity distributions are calculated in a Monte Carlo 
Glauber model [27l-[30|. Both U-l-U and Au+Au follow a 
similar trend for v^/ei and a turn-over is observed in cen¬ 
tral collisions. The inset shows the same quantity but with 
the eccentricity calculated in a constituent quark Glauber 
model [ 1 ^ with the Woods-Saxon parameters proposed 
in Ref. [30|. 


in Ref. (30| with a smaller diffuseness and smaller de¬ 
formation parameter /?2 in combination with the same 
Glauber model, reduces the downturn in central U-l-U 
collisions somewhat but introduces a mismatch between 
the U-l-U and Au-|-Au curves with the Au-|-Au curves 
higher while f2{4}/£2{4} for U-l-U still exhibits a down¬ 
turn (not shown). In the inset of the figure, we show 
the result for a new Glauber calculation using contituent 
quarks as participants 18, JH and the new set of param¬ 
eters [ 3 ^ . This estimate for £2 leads to a seemingly more 
natural behavior for V 2 /S 2 with the drop in the highest 
multiplicity collisions almost entirely gone. The model 
will be investigated and discussed further below. 

The trends of V 2 versus dNch/drj are mostly dominated 
by the elliptic shape of the overlap region in collisions 
with a non-zero impact parameter. To study body-body 
or tip-tip collisions we investigate nearly fully overlap¬ 
ping collisions with minimal activity in the ZDGs. If 
body-body collisions produce smaller multiplicities than 
tip-tip collisions, we expect to see a negative slope in 
V 2 vs multiplicity for these collisions. A negative slope, 
however, can also come from contamination from larger 
impact parameter collisions. To assess their contribu¬ 
tion we use collisions of more spherical Au nuclei as a 
control sample. Figure [3] shows the elliptic flow W 2 { 2 } 
of all charged particles as a function of the normalized 
multiplicity (Mult/(Mult)) for two different systems. We 
increase the acceptance to I 77 I < 1.0 to reduce multiplicity 
fluctuations. The upper panel shows the results for the 
1 % most central events based on the smallest signal seen 
in the ZDGs. Both Au-|-Au and U-l-U show a negative 
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FIG. 3. (Color online) Top panels: charged particle n2{2} vs. 
normalized multiplicity within \r]\ < 1.0. The upper panel is 
for the top 1% most central events based on the smallness of 
the ZDC signal, while the middle panel is for the top 0.125%. 
Small boxes indicate the possible range of variation of V 2 from 
uncertainties in the efficiency corrections on the x-axis. Model 
comparisons are described in the text. Bottom panel: The 
slopes as a function of increasingly tighter ZDC centrality 
selections. The systematic uncertainties are shown as bands. 


slope, which indicates the effect of the impact parame¬ 
ter is still prominent (otherwise we expect the Au-|-Au 
slope to be nearly flat or even positive). The middle 
panel of Fig. [3] shows the 0.125% most central events. 
The negative slope for Au-t-Au collisions is smaller in 
magnitude, indicating the effects from non-central col¬ 
lisions are reduced and the variation in multiplicity in 
Au-|-Au collisions is mainly driven by fluctuations. The 
bottom panel of Fig. |3] shows how the slopes extracted 
from V 2 vs normalized multiplicity evolve with succes¬ 
sively tighter ZDC sections. While the slope for Au-|-Au 
collisions becomes less negative, the slope for U-l-U colli¬ 
sions becomes steeper as the centrality selection is tight¬ 
ened. This demonstrates that the variation of multiplic¬ 
ity in the 0.125% U-l-U collisions is dominated by the 
different geometries made possible by the prolate shape 
of the uranium nucleus and that tip-tip collisions produce 
more multiplicity than body-body collisions. Systematic 
uncertainties shown as bands on the slope were estimated 


by varying the fit range and efficiency corrections. Other 
sources of systematic error are smaller and sub-dominant 
compared to the variation due to the range of efficiencies 
used in the error analysis. Due to large statistical errors, 
no conclusions could be drawn from studies of U2{4} ver¬ 
sus multiplicity in these events. We also measured U3{2} 
in central collisions and found that U3{2} in the 0.125% 
most central collisions are (1.410±0.006) x 10“^ for U-l-U 
and (1.380 ± 0.008) x 10“^ in Au-|-Au collisions (statisti¬ 
cal errors only). The slope of V 3 vs multiplicity was small 
and negative in both systems at about —0.005 ± 0.002. 

The U-l-U data in the top panels of Fig. |3] are com¬ 
pared to the Glauber-cchard model (asssuming V 2 = 
£ 2 (^ 2 )/{£ 2 ))■ The ZDC response was modeled by calcu¬ 
lating the number of spectator neutrons from the Glauber 
model (accounting for the charge to mass ratio of the 
nucleus) and folding each neutron with the known ZDC 
resolution for a single neutron. The Glauber-Xhard model 
significantly over-predicts the observed slope for U-bU. 
This indicates that the variation in multiplicity between 
tip-tip collisions and body-body collisions is smaller than 
anticipated if multiplicity has a significant contribution 
proportional to A^bin- Given this failure, we investigate 
two alternatives with no explicit Albm dependence: a 
constituent-quark Glauber model (Glauber-CQ) [mu 
and the IP-Glasma model [l7| based on gluon satura¬ 
tion [ 1 ^. The Glauber-CQ model neglects Albin and 
counts the number of participating constitutent quarks 
Acq with each nucleon being treated as three constituent 
quarks distributed according to p = poexp(—ar) with 
a = 4.27 fm—1 [l^. This model with = 9.36 mb pro¬ 
vides a good description of transverse energy and multi¬ 
plicity distributions at RHIC [l^ and a better descrip¬ 
tion of V 2 fluctuations than a nucleon based Glauber 
model [ 2 ^. In our simulation, for each Nqq, we sample 
an NBD with parameters tuned to match the distribu¬ 
tions from p-bp 25| and Au-bAu at 200 GeV {n = 0.76, 
and k = 0.34 for \ri\ < 0.5 and n = 2.9 and k = 0.86 for 
1 77 1 < 1). For both Glauber models we use two sets of pa¬ 
rameters for the nuclear geometry, one corresponding to 
the more commonly used values [29| (dashed lines) and 
the new parameters proposed in Ref. (solid lines). 
The effect of the different parameter sets is small. The 
IP-Glasma and Glauber-CQ model are also compared to 
the Au-bAu data (Glauber-ihard is left off for clarity) but 
because of significant uncertainty in the actual shape of 
a Au nucleus, it is difficult to draw conclusions from this 
comparison. 

In U-bU collisions, both the IP-Glasma model and the 
Glauber-CQ model predict slopes closer to the data. In 
the Glauber-CQ model, even though there is no depen¬ 
dence on A^bin, the average number of quarks struck in 
a nucleon (Acq/A^part) is larger in tip-tip than in body- 
body collisions so that tip-tip collisions create more mul¬ 
tiplicity. This leads to a strong anti-correlation between 
A^cq/A^part and £2 which in turn translates into a nega- 











6 


tive slope in V 2 vs. multiplicity. The IP-Glasma model 
exhibits similar behavior. In gluon saturation models 
like the IP-Glasma model, the multiplicity depends on 
QlSj\_/as{Q^s) [l7| where Ql (the saturation scale) is de¬ 
termined by the thickness of the nucleus along the beam 
axis, S±_ is the transverse size of the overlap region, and 
as is the strong coupling constant. For tip-tip collisions, 
the increase in in the numerator will be balanced by 
a decrease of S_\_. In the denominator, however, as de¬ 
creases logarithmically with leading to an increased 
multiplicity in tip-tip collisions compared to body-body 
collisions. 

The slope of V 2 vs. multiplicity provides a detailed 
probe of the multiplicity production mechanism and the 
degree of coherence in nuclear collisions. We find that 
accounting for the observed slope seems to require mod¬ 
els that include effects from sub-nucleonic structure and 
significantly more coherence than is expected from the 
Glauber-Xhard model. Previous studies questioned the 
relevance of iVbin because of the apparent lack of an en¬ 
ergy dependence to Xhard and because the Glauber-GQ 
model also provides a good description of multiplicity 
data. This study however, provides direct evidence con¬ 
tradicting the Glauber-Xhard model. 

In summary, we measured V2{2} and U2{4} for mini¬ 
mum bias, and nearly fully overlapping Au-|-Au and U-l-U 
collisions at yJsjqjq = 200 and 193 GeV respectively. The 
knee structure in high multiplicity U-|-U collisions pre¬ 
dicted by a Glauber model with a two component multi¬ 
plicity model with a dependence on A^bin is not observed 
in V 2 versus dNch/drj. Also, V 2 scaled by £2 from this 
model is found to saturate and then decrease for the most 
central U-l-U collisions. These findings indicate a weak¬ 
ness in the two-component multiplicity calculation that 
is commonly used as part of Glauber models in heavy 
ion collisions. We also used the STAR ZDGs to select 
nearly fully overlapping collisions and showed that for a 
stringent 0.125% ZDG selection criterion, the variation 
of V 2 with multiplicity in U-l-U collisions is dominated by 
the different geometries arising from the prolate shape of 
the uranium nucleus. This demonstrates that ZDGs and 
multiplicity can be used to select tip-tip or body-body 
enriched event samples. The variation of V 2 with mul¬ 
tiplicity in nearly fully overlapping collisions was shown 
to again disfavor the Glauber model including a frac¬ 
tional contribution of A^bin to multiplicity. Models with 
no explicit A^bin dependence such as a gluon saturation 
based model (IP-Glasma) or a constituent quark Glauber 
model agree better with the data. In addition to reveal¬ 
ing fundamental information about the nature of parti¬ 
cle production in heavy-ion collisions, the findings in this 
letter lay the groundwork for more extensive studies of 
the effect of the initial geometry on other observables in 
nearly fully overlapping collisions. 
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